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INTERACTION OP A FREE FLAME FRONT WITH A TURBULENCE FIELD 1 

By Maurice Tucker 


SUMMARY 

Small-perturbation spectral-analysis techniques are used to 
obtain the root-mean-square flame-generated turbulence veloci- 
ties and the attenuating pressure fluctuations stemming from 
interaction of a constant-pressure flame front with a field of 
isotropic turbulence in the absence of turbulence decay processes. 

The anisotropic flame-generated turbulence velocities are 
found to be of about the same intensity as those of the incident 
isotropic turbulence, the lateral turbulence velocities being always 
lower, but the longitudinal velocity is somewhat increased for 
flame-temperature ratios over 7. The small-perturbation anal- 
ysis indicates that the incremental turbulent flame speed is a 
second-order quantity composed of two parts. One part repre- 
sents the root-mean-square area of the turbulent flame front; 
the other represents the contribution of the transverse velocity 
fluctuations resulting from the flame-front distortion. Directly 
at the flame front, the noise-pressure levels of the pressure 
fluctuations are fairly intense (59 to 81 db referred to 0.0002 
dyne/sq cm) even at moderate approach-flow turbulence inten- 
sities. 

INTRODUCTION 

Development of high-output jet engines has stimulated 
interest in the role played by turbulence in combustion 
phenomena. In the earliest studies of flame-turbulence 
interaction, Damkohler (ref. 1) and Shelkin (ref. 2) utilized 
mixing-length theories of turbulence to obtain semiquanti- 
tative relations for predicting flame speeds. Damkohler 
introduced the concept that turbulence of a scale large 
relative to the flame-front thickness increases the average 
flame speed by increasing the instantaneous flame surface 
area. The relations of references 1 and 2 were not confirmed 
by the experimental values of turbulent flame speed in 
Bunsen burners obtained in reference 3. Experiments on 
flames stabilized in channels (ref. 4) suggested that approach- 
flow turbulence had little effect on burning velocity and that 
the disturbances affecting turbulent flame speed were 
primarily flame-generated. A similar conclusion was drawn 
in reference 5. 

In an attempt to obtain agreement between theory and 
experiment, Karlovitz, Denniston, and Wells (ref. 6) and 
Scurlock and Grover (ref. 7) have incorporated the concept 
of flame-generated disturbances in their recent theories of 

* Supersedes NAOA TN 3407, “Interaction of a Free Flame Front with a Turbulence Field, 


turbulent flame speed which utilize G. I. Taylor’s one- 
dimensional theory of diffusion by continuous movements. 
The somewhat arbitrary assumption is made in the analyses 
of both references 6 and 7 that these flame-generated dis- 
turbances constitute additional turbulence only. In refer- 
ence 6 the energy of the flame-generated turbulence is taken 
as the difference between the kinetic energy of the burned 
gas in the absence of turbulence and the kinetic energy 
obtained by using the average velocity of the burned gas 
normal to the turbulent flame front. In reference 7 the 
flame-generated turbulence energy is obtained from a 
momentum balance of unburned and burned gases before and 
after an assumed mixing of the burned gas. 

The data obtained in reference 8 on pentane-air flames 
baffle-stabilized in a rectangular duct suggest that the 
methods of references 6 and 7 considerably overestimate the 
turbulence generated by flame-turbulence interaction. Apart 
from the question of validity of such methods of calculating 
flame-generated turbulence, objections have been raised 
(ref. 9) to calculations of flame speed made on the basis of a 
hypothetical upstream turbulence compounded from 
approach-stream turbulence and turbulence generated down- 
stream of the flame as was done in reference 7 and implied 
in reference 6. 

The present analysis is primarily concerned with the 
turbulence velocities and other fluctuation quantities associ- 
ated with the linearized interaction of a free flame (not influ- 
enced by boun ding walls) with turbulence present in the 
combustible mixture. Such turbulence will be referred to as 
approach-flow turbulence. The flame is treated as a dis- 
continuity specified by the appropriate fundamental (lam- 
inar) flame speed and flame temperature. The interaction 
of such a flame front with a transverse plane wave, that is, 
a vorticity wave or shear wave, of arbitrary inclination rela- 
tive to the front is first analyzed. The effects of an entire 
spectrum of transverse plane waves constituting a weak field 
of turbulence are then developed from the single-wave 
results. The statistical or root-mean-square fluctuation 
quantities describing the pressure fields and the anisotropic 
flame-generated turbulence resulting from interaction of the 
flame front with isotropic approach-flow turbulence are 
obtained for the limiting case of constant-pressure combus- 
tion. Some discussions of turbulent flame speed and of 
combustion noise are also presented. 

” by Maurice Tucker, 1955. 
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FLAME— TURBULENCE INTERACTION PROCESS 

Turbulent motion may be regarded as a Fourier super- 
position of a very large number of different-sized and 
randomly oriented component, plane-wave motions. The 
customary assumptions concerning the turbulence (see ref- 
10) are also made in the present analysis, namely, that 
turbulent decay effects are negligible and that the density 
fluctuations associated with the turbulence are also small 
enough to be neglected. The first assumption, which 
implies inviscid flow and very small turbulent velocity 
fluctuations, permits linear superposition of the component 
waves. With the second assumption, the continuity equa- 
tion requires that these Fourier waves be transverse plane 
waves, that is, vorticity or shear waves. For each of these 
waves the local velocity vector A is perpendicular to the 
vector k, normal to the wave front. The vector k is termed 
the wave-number vector; its magnitude k is termed the wave 
number, which is defined as 2x divided by the wavelength. 
All symbols are defined in appendix A. Any one of the 
parallel planes containing both the local velocity vector A 
and the wave-number vector k is called the “polarization 
plane.” 

Because of the assumed linear superposition of the com- 
ponent waves, the complete interaction results can be ob- 
tained from the study of the interaction of a plane flame 
front with a single-component transverse wave of the turbu- 
lence field. For simplicity, this typical vorticity wave will 
first be taken as a two-dimensional wave. Generalization 
to the three-dimensional case will be made later. The con- 
figuration considered is shown in sketch (a) : 



A flame front moves with mean velocity XJ T into an inviscid 
combustible mixture. This mix ture is at rest, but contains 
a vorticity wave with velocity vector A (components v^, c 
and v',, e ) and with wave-number vector k inclined at an angle 
<p to the positive direction of the a^-axis. In the absence of 
any perturbation interaction, the plane front propagates into 
the combustible mixture with velocity U (the laminar flame 
speed). As indicated in appendix B, the motion of the 
burned gas, whose velocity V is constant, is away from the 


flame front. The flame front is assumed to be of infinite 
extent in directions transverse to the direction of the a^-axis. 

A weak inviscid disturbance field may be resolved into a 
stationary component and a" moving component, both rela- 
tive to the mean local flow (refs. 11 and 12). The moving 
component is an irrotational isentropic pressure-velocity dis- 
turbance. The stationary component, which is conveetod 
by the mean local flow, is a constant-pressure disturbance 
containing any vorticity fluctuations and entropy fluctua- 
tions present in the disturbance field. Thus, the interaction 
of the flame front with the vorticity wave would be expected 
to generate both an irrotational disturbance with velocity 
components u' pyi , v' Pyb and a rotational disturbance with 
velocity components u,, b , v,, b in the burned gas, and an 
irrotational disturbance with velocity components u' P , c , v' p ,c 
in the combustible mixture. The resulting velocity fluctua- 
tions, which include both the irrotational and rotational dis- 
turbances, are designated as u cy v' c and u' b , v' b for the combus- 
tible mixture and the burned gas, respectively. The flame 
front is displaced by an amount £' (y,t) from its mean or 
unperturbed position as a result of the interaction. 

The diagrams of figure 1, which are similar to those used in 
reference 12, may prove helpful in visualizing the interaction 
process. Suppose that at some instant h the flame intersects 
a front of the vorticity wave at point Pi of figure 1 (a) . At 
a later time h+H, the flame has moved a distance U&t and 
now intersects the stationary vorticity wave front at point 
P 2 . A vorticity wave with front parallel to line QP 2 is then 
produced in the burned gas. A cylindrical sound wave is 
generated at point Pi at time U and propagates at speed a b 
into the burned gas while being convected with velocity — V. 
Another cylindrical wave is generated at point Pi at time 
ti and propagates into the combustible mixture -with speed 
a c . The cylindrical wave fronts thus generated form enve- 
lopes (Mach lines) in both the combustible mixture and the 
burned gas, which constitute plane sound waves. 

For the wave-inclination angle <p shown in figure 1(b), 
an envelope cannot be formed on tbo burned-gas side of 
the flame front. The cylindrical sound waves thus 
expand independently and are thereby attenuated. On the 
combustible-mixture side of the front, the cylindrical sound 
waves meet at the common tangent point P 2 . For inclination 
angles less than the critical angle shown, attenuating pressure 
waves are also produced in the combustible mixture until 
another critical angle 180°—^ is reached. Below this 
second critical angle, plane sound waves are again obtained. 
These critical angles may be obtained from the geometry of 

figure 1(b). For small flame Mach numbers where Ms— i 

^o- >i: =sin _1 M 


Vcr, &= sin 
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As the flame Mach number M decreases, attenuating waves 
are produced for a wider range of inclination angle. In the 
lim iting case of very slow flow (constant-pressure combus- 
tion), 3 only attenuating pressure waves appear in combina- 
tion with the vorticity waves if 0°<(<p-O80 o . Quantities 
associated with the pressure wave vanish when <p= 0° or 180° 
because the incident vorticity wave then passes through the 
combustion front without distorting the front. 



(a) Interaction for wave-number vector inclinations generating plane 

sound waves. 

(b) Interaction for wave-number vector inclinations generating non- 

coalescing cylindrical sound waves. 

Figure 1 . — Wave formation arising from interaction of flame front 
with vorticity wave. 

SINGLE- WAVE ANALYSIS 

TWO-DIMENSIONAL FORMULATION 

Tho interaction process described in the preceding section 
is now formulated analytically for the passage of a combus- 
tion front through a single weak two-dimensional vorticity 

• It can bo shown that the static-pressure ratio across a flame front fa given by 

.... 


wave of constant density inclined to the flame front. The 
case of a vorticity wave in three dimensions is considered 
later. 

The combustion front is assumed to be completely speci- 
fied by its laminar flame velocity U and the ratio of stagna- 
tion temperatures r in the burned gas and in the combustible 
mixture, respectively. In the absence of any perturbations, 
the equations for conservation of momentum, energy, and 
mass-flow rate, respectively, as written for a reference frame 
moving with the flame front, are 


Pc+PcU^fr+PtQU- Vf 
y c+! ^ 2 )=~f T „+ 1 (I/-F) 3 
P U= Pb {U-\r) 


( 1 ) 


where subscripts c and b designate stations in the combustible 
mixture and in the burned gas, respectively. For simplicity 
it has been assumed that the combustion process does not 
entail a change in the number of moles per unit mass of gas; 
also, differences in the ratio of specific heats for the burned 
and unburned gas are ignored. The quantity (r— 1) is then 
indicative of the increaso in stagnation enthalpy or heat 
release. 

For the interaction problem the resulting flame-front 
distortion £'(y,t) must be considered in addition to the 
generated disturbances previously mentioned. Thus, both 

the flame-front perturbation velocity and the instan- 

ce' 

taneous flame-front slope will appear in the equations 

of motion. The conservation equations may still be applied 
in a coordinate system moving instantaneously with the 
distorted flame front since extreme gradients occur across 
the front and small disturbances are postulated. The 
various perturbation quantities (designated by primes) are 
assumed to have zero space or time averages. The flame 
speed U T (see sketch (b)) will thus include any time- 



independent contributions arising from the perturbations. 
Conservation of normal and tangential momentum, energy, 
and mass-flow rate provides the following relations: 
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(^+^)[i+(O 2 ]+(Pc+p0[07 I '+^-^) a +(^) a (O , +207r+^-<)i>:s;] 
=to+i>i)[i+(rj a ]+(p 6 +p;)[(^r-Fr+^-<) a +(^) a (f;) a +2a7 r --F 7 .+$;-w;)®u;] 12a) 

(U r +£« — lOf* — ®e=(f7 r — Fr+£«— ' v'b (2b) 

(r+rO (T c +r c )+l QJ T +£-u' c y+\ (pD 2 ]=^- (T b +r b )+\ (u T -v T +t’-u b y+\ (»;)* (2c) 

(Pc+Pc)(Lr+f« y'c-\-v' c ^' p }={pt,+ T —V T+?» — Wi + PUy) (2d) 


Small-perturbation techniques are used to make the interaction problem amenable to analysis. Then, if the velocity 
perturbations are assumed to be small relative to the flame speed U T and the flame-front slope £(, is also assumed to be 
very small so that terms of the second order in the fluctuation quantities may be neglected, application of the linear- 
v' o' T ' 

ized state equation S-= p _ ... | - — and utilization of equations (1) permit the following boundary conditions at the flamo 


P p T 

front to be obtained from equations (2): 


where 


£ +B t ^+B, 
V» 


pt> 


® b n fe i D (i 'U'c 
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U = U +JL *' 
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(3a) 

(3b) 

(3o) 

(3d) 

(4) 


It has also been assumed that the flow upstream of the 
flame is isentropic. It will be shown in the section TUR- 
BULENT FLAME SPEED that TJ t = U is correct through 
first-order terms. 

Another relation is required at the flame front. For the 
two-dimensional case under consideration, the local instan- 
taneous normal propagation velocity U+5U of the dis- 
torted flame front into the combustible mixture at rest (see 
sketch (b)) is 

TTIcTT ZJ T +^i~Uc-\-v' c ^ t 

+ VT+&? 

The incremental propagation velocity SU will be determined 
from existing information on laminar flames. Some sup- 
port for this procedure is given in reference 13. As reported 
therein, radiant flux-intensity measurements on laminar and 
turbulent propane-air flames suggest that a small surface 
element of a turbulent flame is chemically and physically 
the same as that of a corresponding laminar flame. 

The propagation speed of a laminar flame is affected by 
both the ambient pressure and temperature. Although 
the functional relations have not been rigorously determined, 
preliminary indications are that the pressure effect is much 
smaller than the temperature effect. In the present analysis 


the local flame-propagation speed is assumed to depend only 
upon the temperature of the combustible mixture, that is, 



With the empirical relation obtained in reference 14 as a 
guide, it is assumed that U=r 1 +r s 3 r ? where r lf r 2 , and n are 
constants which depend upon the fuel and oxidant under 
consideration. Thus, 

(w) : 

and 

U ^l+A i a =L r fl-^ — u£+ v' e £i 

Correct through first-order terms wherein U T =U, tho fol- 
lowing boundary condition is obtained at the flamo front: 

f, < 3 °> 

In a coordinate system fixed in space, the equations of 
motion for the two-dimensional fluctuation quantities in 
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the burned gas with terms of the second order neglected are 


dt +V 

dpi , T T dp'b 
dt ' dx c 


PbC c 


dn_ 

dt 


-pt,ciV 


7 dlll_ 

1 ^ 


dx e ' 

Pb dx e 


rdv't,_ 

IM 


dx e 

pb dy 


Apt, ( 

'duj dvj'y 
\dxc dy y 

)=0 

■ dT' b 
dx e 

. (dUb 

f Hs; 4 

dvi\ 
d y) 


For a coordinate system moving with constant velocity V, 
the preceding equations reduce to the same form as the 
corresponding equations for the fluctuation quantities in 
the combustible mixture relative to a coordinate system 
fixed in space. Thus, the flow equations for both the 
combustible-mixture fluctuations and the burned-gas fluc- 
tuations may be written, with appropriate subscripts c or 
b, as 


du' 1 dp^ 

dt~ p dx 



(5) 


where v—k[x 0 J r k' i y, and k[ and k'% are components of the 
wave-number vector k in the x 0 - and y-directions, respec- 
tively, with ki/k[=tB,n <p. As a result of the linear boundary 
conditions of equations (3), obtaining a unique solution of 
the interaction problem requires that the arguments of all 
disturbance waves match at the flame front. This matching 
requirement together with provision for differences in phase 
angle yields the following form for the vorticity wave 
present in the burned gas: 

where 

<k - Kxb+fyj 

Pressure fluctuations generated by the interaction must 
satisfy the following wave equation, with appropriate 
subscripts c or b, which is obtainable from equations (5) : 



The present analysis will be concerned with the limiting case 
of very slow flow (constant-pressure combustion). It is 
clear from the relations given for the critical wave-inclination 
angles c and <p a ,b that, for very slow flows, only the irrota- 
tional isentropic pressure waves described in the section 
{TAME-TimBULENCE INTERACTION PROCESS will 
be generated by the interaction. The form of these pressure 
waves that satisfies equation (6) has already been established 
in reference 12 in terms of the variables y and f , where 


The coordinate systems x c ,y and x t ,y for the first-order 
analysis are indicated in sketch (c) . 


/* o /=■/ 



A two-dimensional vorticity wave in the combustible 
mixture with velocity vector of magnitude A that has its 
wave-number vector k inclined at an angle <p to the positive 
direction of the a^-axis may be written in the form 


nc^(Jc)(x-Ut) 

bcXc+Ccy+Udct 

Vb=ifb)[-x b +(U-V)t] 

fi, == hi,x b -\-Ct,y-\-U dbt „ 


The variable y is proportional to the distance from the flame 
front. At the front, i 7 „=i?ft= 0 ; upstream of the flame, 
y c is positive; and downstream of the flame, is positive. 
The equation constant defines planes moving with con- 
stant velocity ( U<Z)c,b at an angle tan -1 ( c/b to the flame 
velocity U. Equation (6) takes the form of the Laplace 
equation 



when 


/dyv /d£v = 2 

a LW \dz/ \dy) J 

dy d£ . dy 5£ 
dt dt dx dx 


( 8 ) 


<c 

U 


n 

{A sin <p)e ir , -jj=(—A cos <p)e u 


Matching arguments of the pressure and vorticity waves at 
the combustion front where y=0 and satisfying the require- 
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ments of equations (8) provide the following values for the 
constants of equations ( 7 ): 


k[M- , 


(U-V'P TP ’ 

\ TJ ) al 


4 = 




1 -M s 


d b =- 


K 


bi+p-dur- 


/Ll M | 

( u-vyzp’ Jc ~ i—M 2 

\ U ) al 




TP 

bt+ct-dl hL 

at, 


( 9 ) 


To satisfy the requirements that the arguments of all dis- 
turbance quantities match at the flame front and that 
d*f d 3 £ 

— -^=-— the flame-front slope must take the form 
otoy byot 

^=(B r 1 +iH J )( tan <p)e u (10k) 

Integration of equation (lOj) with respect to time t gives the 
flame displacement as 

r=^7 {irP-iEpy* (io t) 

h, j /tj 

EXTENSION TO THREE-DIMENSIONAL DISTURBANCES 


In addition to the boundary condition from equations ( 3 ), 
the pressure fluctuations will be required to satisfy the 

boundary condition — =0 at 17= <° . Utilizing equations 


( 36 ) and ( 37 ) of reference 12 yields 


£=(R 1 +iR 2 )e , *‘~’'° 

Pc 

(10a) 

2‘=(J 1 +U 2 )^-^ 

pb 

(10b) 


In the combustible mixture, density fluctuations are associ- 
ated only with the pressure fluctuations according to the 
/ 1 w 

isentropic relation — =- — • In the burned gas, density 

Pc ‘ Pc 

fluctuations may also be caused by entropy fluctuations 
generated by flame distortion and heat-release fluctuations, 


if present, as well as by pressure fluctuations. Velocity 
fluctuations are associated with both pressure fluctuations 
and vorticity fluctuations. It is convenient to deal with 
the pressure coefficient p'/p c U 2 . Thus, the disturbances 
arising from the interaction of the flame front and incident 
vorticity wave take the following forms: 

iRi+iB^—tRO’+W*)*' 

(10c) 

= ( J (1 > + ( 10 d) 

PcU PcU 2 

-=4 {Ji+iJJe^+iLt+iLtW* 
Pb 1 

(10e) 

^(Nt+iNJe^+iGt+iGde 1 * 

(10f) 

^=(P 1 +tP 2 y r ‘-’‘+(/>+ ij 2 y+ 

(10g) 

^=(n ? i+ iW^e iic ~’’ c -\-{A sin <p)e u 

(10h) 

jj=(X 1 +iX 2 )e iCc ~' ,c —(A cos <p)e u 

(lOi) 

The flame displacement velocity may be written 


^=( 27 , + = (Hr+iH.y* 

(10j) 


Equations ( 10 ) describe the interaction of a flame front 
with a constant-density vorticity wave having velocity com- 
ponents 


sin <p)e ir and v -jj—{—A cos tp)e t ’’ 

in the x c - and ^-directions, respectively. The vorticity wavo 

may be considered to have a third velocity component 
/ 

-jj=Ce u in the z-direction. In the preceding linearized 

analysis, the amplitude G was not prescribed. This com- 
ponent, which is normal to both the u’ and v’ components 
and parallel to the plane of the flame front, then is associated 
with a corresponding component of the vorticity wave in the 

burned gas ^j?=Ce { *. 

Inasmuch as turbulence fields are three-dimensional, the 
interaction equations must be revised accordingly for applica- 
tion to the spectral analysis which follows. Assume, as 
shown in sketch (d), that the polarization plane which -will 
contain the wave-number vector k is inclined at some angle 6 
to the z^jZj^-plane of a new coordinate system x Ut , x 2 , e , x 3 , e 
fixed in space for the combustible mixture. 



The corresponding coordinate system x Ui , x 2>b , x 3 , b for the 
burned gas is assumed to be moving -with mean flow velocity 
V. Components k u k 2 , and k t of the wave-number vector k 
in the directions of x hc , x*, c , and x 3 , c , respectively, are 


ki=k cos <p 
k 2 =k sin <p cos 6 - 
k z =k sin <p sin 8 „ 


(ID 


With primed perturbation vector quantities referring to the 
original coordinates x c , x b , y, z and unprimed perturbation 
vector quantities referring to the x 1 , t , x 2 , c , x z , c - and x ub , x 2 , b , 
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Xz, s-coordinate systems, the following transformation rela- 
tions apply: 


uue=u e 

cos 8 — w' c sin 8 'ih, b =v' tl cos 8 — w' b sin 8 

^xc—v'e, sin 8 -\-w' c cos 8 %6=vs sin 6 - {-w' b cos 9 


f=r 

cos 6 

h u =5i sil >« 


y (i 2 ) 


This notation refers only to equations (12) . 

With the use of equations (10) and (12), the interaction 
fluctuation quantities (again designated by primes) referred 
to the coordinate axes x UC) a^, c , Xz, c and x 1)b , Xt, b , X 3 , b may be 
written as 


J^ = (RW + iRQ)) e <te-*c 


(R.+iR^e 1 ^ 

Pc 7 

^=- (JL+iJJe'^+ih+ihy* 

Pit 7 

( Wj + iWi) e ,rc_ ’ ,c -|- (A sin <p)e u 

t 

l ~=(X l J r iX a) (cos 0) (A cos ^cosO+C sin 6)e ir 


q j~=(Xi+iX 2 )(sia 8 )e , ^ c ~ Vc —(A cos tpsind—G cos 8 )e ir 


r (13) 


^A=(Ni + iNi) ((?i + iGj) e** 


solution for the attenuating pressure-wave regime is indicated 

in appendix B. Inasmuch as the flame Mach number M^~ 

d C 

is generally much less than 0.01, the limiting case of very 
slow flow (constant-pressure combustion) provides a reason- 
able simplification of the problem and only the attenuating- 
wave solution need be considered. The amplitude coeffi- 
cients for this limiting case are given in appendix C (eqs. (C8)) . 

Equations (13) and (C8) describe the linearized interaction 
of a constant-pressure flame front with a single vorticity 
wave or shear wave having its wave-number vector k inclined 
at an angle <p to the direction of travel of the undisturbed 
flame front and having its plane of polarization inclined at 
an angle 8 to the x Uo , a^, c -plane of the coordinate axes x 1>c , 
Xt ,c, Xi, c . The wave-number vector of the shear wave 
generated in the burned gas makes an angle p 6 =tan -1 
(r tan ip) with the direction of propagation of the undisturbed 
flame front. Attenuating potential fields are generated in 
both the combustible mixture and the burned gas. Physical 
quantities associated with these fields attenuate exponentially 
with increasing distance from the flame front. The ampli- 
tude coefficients for a given heat release (a prescribed r) 
and a given inclination angle <p and a polarization angle 9 
depend upon both the intensity of the incident vorticity 
wave and the heat-release perturbation parameter t'/t. In 
the absence of such heat-release perturbations, there are no 
density fluctuations in the burned gas (correct to order AF). 

These single-wave results may be used to determine the 
interaction of a constant-pressure combustion front with a 
turbulence field of constant density for the case of negligible 
turbulence decay. The turbulence field will contain an in- 
finite number of transverse plane waves with all wavelengths 
and planes of polarization. The spectral analysis technique 
used in obtaining such a superposition of waves will be dis- 
cussed briefly before proceeding with the interaction problem. 

SPECTRAL ANALYSIS 

GENERAL CONSIDERATIONS 

A turbulence field satisfying the incompressible-flow con- 
tinuity equation may be represented by the following super- 
position of plane transverse waves: 


■^?^=(Pi+iPj)(cos0)e f6 ’* + [(A+ i A) cos 0 — O sin 0] ^ 
^^==(Pi+iPj) (sin 0)«5 ,r6_,6 -(-[(/i+'tJa) sin 8-\-G cos 8} e** 

~={H l +iH J )e u 

(Hj+tSaX tan <p cos 8 )e u 
c =(Z7 1 +'iiJ 2 )(tan <p sin 8 )e <a 

5! J 

The fluctuation amplitude coefficients of equation (10) or 
(13) may be determined from equations (3) and (5). Details 
of the solution are given in appendixes B and C. The general 


JIf «- 2 dZ(t,f) 

OD 

where x is a position vector, k is a wave-number vector, t is 
the time, and d ZQc,t) is the random amplitude vector of a 
component wave. The quantity d Z(k,t)eAl represents the 
contribution to the velocity field from a volume element dt 
in wave-number space. When, as in the present case, the 
equations of motion are linear there is no modulation or 
interference between component waves, and the various sta- 
tistical quantities describing a random field may be obtained 
from the results of a single-wave analysis. To avoid the 
interpretative difficulties associated with the random variable 
ZQc,t), which is nondifferentiable with respect to k, use is 
made of the techniques of references 10 and 15, which utilize 
correlation spectra rather than amplitude spectra in the 
analysis of homogeneous turbulence. 
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A velocity correlation is defined as the ensemble average 
u t (x,t)uj(x' ,t) of the product of a fluctuation-velocity com- 
ponent Ui at x and a component Uj at x' =x-\-r where r is a 
separation vector. The subscripts i and j take on the values 
1, 2, and 3. The ensemble average, designated by a bar, 
may be regarded as the result of averaging the product 
u ( (x,t)uj(x,t) at a given instant over a very large number of 
statistically similar fields. The nine velocity correlations 
uflij constitute the velocity correlation tensor Tu(x,x',t). 
For a homogeneous field, T {) depends only on r so that the 
tensor may be written TyfrJ). 

As shown in references 10 and 15, the velocity correlation 
tensor has the following Fourier integral expansion: 

T (] (r, dJVK&#)“ JJJ (14) 

where F t] (k,f) is the spectral tensor function, <£«(&,*) is the 
spectral tensor density of a homogeneous turbulence field, and 

$,Xk*)d£=dZKp)dZ^) (15) 

where d Z\(k,t) denotes the complex conjugate of dZ t (k,t). 
For r=0 and i—j, equation (14) may be written 

GO 

T«(0, t)=%= JJJ AZT^DdZ^t) (16) 

GO 

For homogeneous turbulence fields, wherein ensemble aver- 
ages and space averages are identical, equations (15) and (16) 
provide the basis for obtaining the spatial mean-square 
velocity components from the single-wave results given by 
equations (13) and (C8). Equations (15) and (16) are also 
applicable to scalar fields. In the absence of viscosity, as 
postulated, the shear-velocity fields present in the com- 
bustible mixture and in the burned gas are homogeneous, 
and application of equation (16) presents no complications. 
The corresponding potential-flow fields, although spatially 
inhomogeneous, are homogeneous in the given Zs,X 3 -planes. 
It has been shown in reference 16 that equation (16), in 
effect, may be applied for such fields to obtain the mean- 
square fluctuations pertaining to a given plane of homo- 
geneity. 

As a result of the preceding discussion, the single-wave 
interaction results for constant-pressure combustion will be 
used to obtain the spectral densities of the fluctuation 
quantities at the flame front where the attenuation factors 
&~”c and e-'a are unity and f 0 =fb=<r=^=p. For concise- 
ness, define BA sin p=A a) -KA (2) ,fFA sin 9 >=Wi-HW 2 , 


GA sin <p=Gi-\-iGi, NA sin #>=IVi+i2Vj, HA sin tp^H-\-iH-i, 
and C n A sin <psC. With the notation 



and, for example, 



as in equation (16), the following equations are obtained by 
analogy with equations (13) and (C8) for the case where 
heat-release perturbations are absent, that is, r' = 0: 

d (^f^WdZ^+dZ,,, 
xW cos0dZi, c — (cot ip cos 6+C 1 sin 0)dZ IfC 
xW sin0 dZ IiC — (cot <p sin 0— C 1 cos 0)dZi iO 
d d Z uc +G dZi.c 

d (j^=xN cos 0 dZ x . c -( c - W T e — Q+ O l sin 0 ) dZ,. , 
d (^f^=xN sin 0 dZi, c -( COt ^ Sm - G-C l cos0)dZ,. e 

d(«=.-fdz,,~ ^dZ,- 

d(|)=HdZ 1>c 

d(£ } 0 )=-ff tan <p cos 0 dZi, e 
d(£* s )=H tan <p sin 0 cLZi,, 

In the velocity ratios of equations (17), the first term on the 
right side represents the potential-flow contributions 

or Ui jj the second term represents the shear-flow contri- 
butions or The subscript i takes on the values 

1, 2, and 3. 
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At a given instant the spatial mean-square potential- and shear-flow contributions of the disturbance fields are obtained 
from equations (17) : 




(^ f )-/// dZf - dZ - 


(^)+0f) 2 = JJJ dZt, e dZ 2 . c +dZf, c dZ, c 


(^)=(^)+(^)=JJJ (W*W) dZf, c dZ lrC 


( ! rr) 2= //J (6r * <?) dZ * • ^ • 


dZf, c d z ue 


(^) 2 = (Tr) 2 + (^r) ,= JJJ 1 

— OO 

The moan-square flame-front quantities are 

© 2= d ^.c dZ,. e 

— CO 
— CO 

CD 

H £+H7= JJJ(27*H) tanV dZ?. e dZ 1>c 


(18) 


From equations (C8) of appendix: (7 with t'=0: 


R * R _ (r l) a (r tanV l ) 2 
A sinV 

tanV— l) 2 
A 

a*n— t 2 secV [A — 4 t(t 2 — 1) tanV] 
u J A(1 +t s tanV) 

r 8 (r— l) 2 (r tanV-l) 2 secV 
A(l+r l tanV) 


E*H-- 


4t 5 secV 


(19) 


MEAN-SQUABE FLUCTUATIONS FOE INITIAL ISOTBOPIC TUEBULENCE 

For a given combustion process (r and Z7 prescribed), the 
spatial mean-square fluctuations of equations (18) depend 


upon the quantity dZf,„ dZ llC , which is specified by the type 
of turbulence present in the combustible mixture. The 
results obtained in reference 16 for the interaction of axi- 
symmetric turbulence with a shock wave suggest that for 
the present problem, the degree of anisotropy of the incident 
turbulence field may not be critical. For simplicity, the 
turbulence in the combustible mixture is assumed isotropic. 

As indicated in reference 10, the spectral density ten- 
sors for any isotropic turbulence field satisfying the 
incompressible-flow continuity equation are 


<iV©=n(*)(^u-^) (20) 


where & 2 =&?+&f+&!; 5 W =1 for i=j ; 5^=0 for and 
Q(k) is the scalar amplitude function defining the spectral 
density tensor. From equations (15), (20), and (11), 


dZf, „ dZj, e = ^(b)^ 2 sin 2 p dk 


43G875 — 57- 


40 
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or, transforming to spherical polar coordinates k, 6, p, wherein 
dkzsdki dk 2 dk z =lr sin p dl* d0 dp, yields 

dZf, c dZ lf c = Q(k)k A d k d0 sin 3 p dp 

dZf c dZ> c = dk d# sin pfcos 2 p+sin 2 p sin 1 0]dp 

V (21), 

dZf, c dZ 3i „= d& d0 sin pfcos 2 p+sin 2 p cos 2 Q]d<p 

dZf , c dZj_ c + dZf , e dZ 3iC = 0 (k)k* d& d0 sin p[l + cos 2 p] dp , 

The mean-square velocity components of the incident isoj 
tropic turbulence field are then given by 

J" Of*)* 1 dk J* T d9 JJ sin 3 p dp=|^” Q(k)k* dk 

(^p)*+ (cu*y = Jo" dt J Q 2r ds J o " sinp(l+cos 2 p)dp ► 

i 

3 Jo J\ 

( 22 ) 

The remaining spatial mean-square fluctuation quantities!, 
referred to the intensities of equations (22) in order that the 
scale need not be specified, are given by 


/«».«. »v i 

V u ) u ) 3 / _ a ^ p T (A+4T8ec a p)sin 3 p 

7x 7A 2 , AwY ACl+AtanV) 

\ru~) + \u~) ( 

Y /W? i 7 %mY 

v a ; UtJjaJLz 

W (V)‘ 


-3 r (r tan 2 P 

4 Jo A(1 


— I) 2 sin p tan 2 p 


A(1+t 2 tan 2 p) 


dp (23e) 


fiiY 

W/ P r sin p tan 2 p , 

erf J " A ' 


m J ° A ' 


(MiM 


kr-iyf^irttmt'p-Wv (23d) (Jtr)* T Jo 


P SlikWdk] ri . , 

Jo r sin 3 p 

P 0/WJ-4JI. Jo A cos 4 p 


Q(Ar)AdA: 


dp (23h) 


aw> 2 ^Y+P^Y 

V U ) \U )^\U J 3, ,„rsin 3 p, , , 

TSUa- 1 KY — < rt “V-i)‘d„ 


A 2 

/« U'.e\* 

) 

\u J 



v u , 

1 3A P T [A 

Ab,.A 

? 4 Jo 

V u , 

) 


A(l + r 2 tah 2 p) 


The subscript 1 designates a longitudinal component; sub- 
2 scripts 2 and 3 designate the lateral components. Of the ro- 
1J dp maining subscripts, s denotes a shear-flow component; p 
(23b) denotes a potential-flow component; c refers to the combusti- 
ble mixture; and b refers to the burned gas. Equations (23a) 
to (23g) have been integrated numerically using Simpson’s 
dp (23c) rule the following increments for tt: 2° intervals from 

v ' ' 0° to 20°, 5° intervals from 20° to 70°, 2° intervals from 70° 

to 90°, and so forth. Numerical results are listed in table I. 


TABLE I.— FLAME— TURBULENCE INTERACTION FLUCTUATION RATIOS 

I 
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DISCUSSION OF RESULTS 

Mean-square fluctuation quantities generated by the 
linearized interaction of a constant-pressure combustion 
front with a weak isotropic turbulence field satisfying the 
continuity equation for incompressible flow are given by 
equations (23) in terms of the incident-turbulence-velocity 
components. Equations (23a), (23b), and (23e) apply only 
at the flame front where the attenuation factors are unity. 
This restriction does not apply to equations (23c) and (23d) 
in the assumed absence of turbulent decay processes. 

VELOCITY FLUCTUATIONS 

Potential-flow fluctuations. — The root-mean-square ve- 
locities associated with the attenuating pressure fields gen- 
erated in the combustible mixture (eq. (23b)) and in the 
burned gas (eq. (23e)) are plotted in figure 2. Since these 



Floras 2. — Effect of flame-temperature ratio on shear-flow and 
potential-flow velocity fluctuations. 

ratios apply only at the flame front where the exponential 
attenuation factors are unity, they represent maximum values. 
As is to be expected from the boundary-condition require- 
ment of equal pressure fluctuations with differing densities 
on each side of the flame front, potential-flow velocity com- 
ponents in the burned gas exceed those in the combustible 
mixture. Both ratios increase with increasing flame- 
temperature ratio, reaching asymptotic values of V3/2 in the 
burned gas and 1 in the combustible mixture. Hot-wire in- 
strumentation will respond to these fluctuation velocities as 
well as to the shear-flow fluctuation velocities. In view of their 
exponential attenuation characteristics, however, such con- 
tributions would not be of importance unless measurements 
were made at stations very close to the flame, that is, within 


a distance of the order of incident-turbulence scale. For 
the low flame Mach numbers encountered in combustion, the 
contribution to the hot-wire signal voltage of the unatten- 
uated sound waves described in the section FLAME — 
TURBULENCE INTERACTION PROCESS should bo 
quite small. 

Flame-generated turbulence. — The shear flow in the 
burned gas (eqs. (23c) and (23d)) constitutes the flame- 
generated turbulence occasioned by the presence of approach- 
flow turbulence. These velocities, referred to the incident 
turbulence velocities, are also plotted in figure 2. A slight 
amplification of the longitudinal component occurs for 
values of r over 7. In the limit, as r becomes very large the 
longitudinal and lateral velocity ratios approach asymptotic 

values of t/ 3/2 and ^ -\fS, respectively. 

The diagrams of figure 1 indicate that a pressure wave 
interacting with the flame front can also bring about a shear 
flow in the burned gas. Although the reflection and the 
consequent impingement of the pressure fields described by 
equation (23a) upon the flame front are possible, any addi- 
tions to the flame-generated turbulence level through the 
reflection process would probably be negligible because of the 
attenuating nature of the pressure field. Thus, contrary to 
the predictions of references 6 and 7 that the flame-generated 
turbulence intensity should be many times greater than the 
intensity of the incident field, the present analysis indicates 
that the two intensities are about equal. 

It is interesting to note that a stream contraction (ref. 17) 
increases the downstream velocity of the mean flow (as does 
the flame front also) while exercising a different selective 
effect upon an incident isotropic turbulence field. For ex- 
ample, with- a sevenfold increase in the downstream velocity 
of the mean flow, the longitudinal velocity ratio (in the ab- 
sence of decay effects) is 1.01 for the flame front and 0.31 
for the contraction. The corresponding lateral velocity 
ratios are 0.87 and 2.29, respectively. 


TURBULENT FLAME SPEED 

The higher mass-flow burning rate of a turbulent flame as 
compared with that of the corresponding laminar flame is 
generally described in terms of a turbulent flame speed U T . 
The flame-speed ratio U T IU is generally assumed to be 
equivalent to the ratio of turbulent-to-laminar flame surface 
area. 

Calculation of the turbulent flame speed requires con- 
sideration of second-order terms. The local instantaneous 
normal propagation velocity 17+517 of the distorted flame 
front into the combustible mixture at rest is 


U+dU- 


U T +^',— ttl.e+fij^4.e+5i i e U3, 






(24) 


Let f7 r = 17+ l/i+ U n -\- ■ ■ ■ where U N represents a steady- 
state contribution to the flame speed of order N. The 
perturbation quantities, for example, are written as 

£=£+£ T + ■ • • 


c — • • • 
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The superscript on a fluctuation quantity indicates the 
order of the perturbation. As before, 5U/U is taken as 


SU nr a Tl T'c (7—1 )nr 2 
U~ U T c ~ 





Substitution of this expression into equation (24) and per- 
forming the indicated expansions yield the following rela- 
tion, which is correct through second-order terms: 


Averaging this equation yields 



U* 

V 


=r 


pJJ* 


'2 [(4,t)'+(S3,«) ! ] V JJ £*S, 0 JJ ^*3, C 


For the limiting case of constant-pressure combustion, r=0. 
Thus, the ratio of turbulent-to-laminar flame speeds can bo 
written for this case as 


r (^M r ^ t( ^* ) ’ +(au ’ 1 } 

' Y ' ' Y 7 

first-order second-order terms 

terms 



first-order 

terms 


£» ^l.C I ‘ c uj I 'llj.C yl \ 

u u 


second-order terms 


^-i+| [(a t .)*+(a l ,)l-(% ! a,„ 


(25) 


The first two terms on the right side of equation (25) 
represent the ratio of averaged turbiflent-to-laminar flame 
surface area. The third term is a correlation coefficient 
representing the transverse-velocity-fluctuation contribu- 
tion to the turbulent flame speed caused by flame-front 
distortion. It is interesting to note that only the trans- 
verse velocity fluctuations appear explicitly. The second 
term has already been determined (eq. (23g)). Tho third 
term may be obtained from equations (C8) as 


<c vi <c _§Z 

U u 2 



[t(t— 1) tan 4 H-OrH-l) tan 2 p+(r+l)] 




(26) 


The imaginary term has not been written out inasmuch as 
it does not contribute to the integral. The results of the 
indicated integrations are listed in table I. 

Equation (25) for the flame-speed ratio U T /U may also 
bo witten in the form 

(27) 

where the flame-speed parameter S is obtained from the 
values listed in table I from equations (23g) and (26). The 
variation of this parameter with the flame-temperature 
ratio r shown in figure 3 suggests, on the assumption that 
the flame-front slopes £r 2 e and govern contributions to 
the right side of equation (25), that the flame front with the 
higher heat release is distorted less by a given intensity of 
turbulence. For the degenerate case r=l (no heat release), 
the parameter S becomes infinite — a condition compatible 
with this viewpoint. 

The present analysis requires that the flame-front slopes, 
as well as the other perturbation quantities, be small. The 
preceding discussion suggests that the incremental flame- 
speed ratio (U T — XT) ITJ may be of the second order as a 
result of this restriction to small flame-front slopes. 



Figure 3. — Effect of flame-temperature ratio on flame speed for 
incident isotropic turbulence. 


COMBUSTION NOISE 

The root-mean-square pressure-fluctuation coefficient 
-V^/pc^ which applies directly at the flame front, is 
plotted in figure 4 for the limiting case of constant-pressure 
combustion. The pressure fluctuations are a measure of the 
random noise generated by the interaction of the flame front 
with the incident turbulence. In acoustical measurements, 
the noise level in decibels is usually given with respect to a 
reference pressure of 0.0002 dyne per square centimeter 


(ref. 18), which corresponds approximately to the pressure 
amplitude of a plane sound wave of minimum audible 
intensity at a frequency of 1000 cycles per second. The 
noise-pressure level in decibels is defined by the relation 


Noise-pressurelevel r 20 logi 0 



=74-f 20 log 10 ^2 ( 28) 


where the pressure fluctuations in the combustible mixture 
are given in dynes per square centimeter. 
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TABLE II.— HYDRO CARBON- AIR FLAME DATA AT SPECIFIED CONDITIONS 


[Statlo pressure and temperature of combustible mixture, 760 mm Hg and 26° 0, respectively.] 



Propane, 0:Hs 

Ethylene, CiHt 

Acetylene, CiHj 

Atstoicblometric 

At mn.Tlmnm XJ 

Atstoicblometric 

At maximum XJ 

Atstoichlometrlo 

At maximum XJ 

Fuel-air ratio 

0.06S8 

0.0721 

0.0676 


A 0753 

a loio 

Fuol In air, porcont by volume 

4.04 

4.54 

6.54 

7.65 

7.75 

10.70 

Laminar fliuno speed, U, cm/sec 

37.5 

39 

64 

es 

123.5 

141 

Adiabatic flame-temperature ratio, t... 

7.70 

7.45 

7.95 

8.00 

8.64 

8.71 



FiauiiB 4. — Effect of flame-temperature ratio on random pressure 
fluctuations generated at flame front. 


Equations (23a) and (28) indicate that the noise level 
should be particularly dependent upon flame speed. Propane- 
air and acetylene-air combustions, which are character- 
ized by a low flame speed and a fairly high flame speed, 
respectively, will be considered for illustrative purposes. 
Pertinent data for these flames at maximum-flame-speed and 
stoichiometric conditions for an ambient temperature of 
25° C and a pressure of 760 millimeters of mercury are given 
in table II. The adiabatic equilibrium flame temperatures, 
at which the total enthalpy of the fuel and oxidant equals 
the total enthalpy of the products of reaction, were calcu- 
lated using the procedure of reference 19. (Total enthalpy 
includes the chemical contributions to the internal energy.) 
Flame-speed data were obtained from references 20 and 21. 

If the flame-front turbulence intensity is assumed 

equal to 10 percent, noise-pressure levels of 59 and 81 decibels 
are obtained for propane-air flames and acetylene-air flames, 
respectively, under conditions for maximum laminar flame 
speed. At an approach-flow velocity of 1225 centimeters per 
second, which is in the range of velocities usually encountered 
in combustion experiments, the corresponding intensity of 
the approach-flow turbulence would be about 0.3 percent 
for the propane-air mixture and about 1 percent for the 
acetylene-air mixture. 


Thus, the pressure fluctuations generated at the flame 
front when the incident turbulence is of low intensity, 
although small compared with ambient pressure, are appar- 
ently of fairly high acoustical intensity for constant-pressure 
combustion. Because of the exponential attenuation of 
these pressure fields, the “far-field” acoustic intensity (the 
intensity at distances very far from the flame front) ap- 
proaches zero. For cases other than constant-pressure 
combustion, a finite “far-field” intensity is obtained. 

CONCLUDING REMARKS 

The present linearized analysis has treated the interaction 
of a field of isotropic turbulence with a free flame front under 
constant-pressure combustion conditions with no turbulence 
decay processes or heat-release fluctuations. The interac- 
tion produces an anisotropic turbulence field in the burned 
gas which has axisymmetry about the mainstream direction. 
Contrary to the results predicted by several current theories 
of turbulent flame speed, the flame-generated turbulence 
velocities caused by approach-flow turbulence do not differ 
greatly from the turbulence velocities of the incident field. 

The incremental flame-speed ratio ( U T — U)/U as obtained 
from the present analysis is a second-order quantity con- 
sisting of two parts. One part represents the root-mean- 
square area of the turbulent flame front; the other represents 
the contribution of the transverse velocity fluctuations which 
result from the flame-front distortion. The flame-speed 
ratio Ut/U for a given level of incident turbulence intensity 

is found to decrease with increasing heat-release 
rates (increasing values of r). 

Random pressure fluctuations generated in both the com- 
bustible mixture and the burned gas, although small com- 
pared with ambient pressure, give rise to appreciable noise 
levels (59 to 81 db) directly at the flame front even for very 
moderate intensities of approach-flow turbulence (flame- 
front turbulence intensities of 10 percent). For the limiting 
case of constant-pressure combustion, the pressure waves 
attenuate exponentially with distance from the flame front, 
so that the “far-field” intensity approaches zero. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 25, 1956 

















612 


REPORT 1277 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX A 


SYMBOLS 


A 


A 


a 

Bi,B a ,B 3 

b 

c 


Cp 

c, 

A 

d 

E 

Fu(k,t) 

f 

G 

<?„<?* 

E 

E u E a 


hl,h -2 . . . hio 

h,It 


Jliji 


K u K a ,K 3 ,E 

k 

k 


Li, Li 

M 
m 2 
N 

N u N 2 


n 


Pl,P 2 


P 

f 

Q 


magnitude of two-dimensional vorticity- 
wave velocity vector in combustible mix- 
ture 

two-dimensional vorticity-wave velocity vec- 
tor in combustible mixture 
speed of sound 
coefficients defined in eqs. ( 4 ) 
constant defined in eqs. ( 7 ) and ( 9 ) 
amplitude of combustible mixture shear- 
wave component parallel to 23, c -axis 
constant defined in eqs. ( 7 ) and ( 9 ) 
specific heat at constant pressure 
specific heat, at constant volume 
WU-V), eqs. ( 4 ) 
constant defined in eqs. ( 7 ) and ( 9 ) 

-TV U, eqs. ( 4 ) 
spectral tensor function 
constant defined in eqs. ( 7 ) and ( 9 ) 
(Qi+iG a )/A sin <p 

amplitude coefficients of shear-wave longi- 
tudinal component in burned gas 
(Hi+iH a )/A sin <p 

amplitude coefficients of flame-front dis- 
placement 

groupings defined in eqs. (B 16 ) 
amplitude coefficients of shear-wave trans- 
verse component in burned gas 
amplitude coefficients of pressure wave in 
burned gas 

J<1) s Sp Jl > jm J » «*• (10d ) 

coefficients defined in eqs. ( 4 ) 
magnitude of wave-number vector k 
wave-number vector of shear wave in com- 
bustible mixture with components k[,k a 
in 2 c ,j/-coordinate system, with com- 
ponents ki, k 2 , k 3 in x uc , x?,e, 23, ,,-coordinate 
system 

amplitude coefficients of density associated 
with shear entropy wave in burned gas 
flame-front Mach number, M=U/a c 
a^/yU 3 , eq. (B 16 ) 

(Ni+iN a )/A sin <p 

amplitude coefficients of longitudinal ve- 
locity component associated with pressure 
wave in burned gas 

exponent used in representation of laminar 
flame speed as function of combustible- 
mixture static temperature 
amplitude coefficients of lateral velocity com- 
ponent associated with pressure wave in 
burned gas 
static pressure 
static-pressure perturbation 
term not contributing to the integral in 
eq. ( 26 ) 


R 

Ri,Ri 

R m ,R m 

T 

ri,r 2 


S 

T 

Pn(z>t) 

T s 

t 

U 

U T 

u' 


iq. cj Uv.ct *^8, e 


1h, b, 'Ih.b, U 3 . b 


V 

V T 

v' 

W 

W u 


w' 


Xi, X2 


x 


x c 




*2.6 


X* .b 


(R m +R<*)IA sin <p 
gas constant 

amphtude coefficients of pressure wave in 
combustible mixture 
R (l >=R,hM\R’v=R a lyM*, eq. ( 10 c) 
separation vector 

constants used in representation of laminar 
flame speed as function of combustible- 
mixture static temperature 
flame-speed parameter, eq. ( 27 ) 
static temperature 

velocity correlation tensor for homogeneous 
turbulence 

stagnation temperature 
time 

laminar or fundamental flame speed 
mean turbulent flame speed 
longitudinal component of velocity perturba- 
tion 

velocity perturbation components in com- 
bustible mixture parallel to x it c ,2o, e ,Xx r 
coordinate axes, respectively 
velocity perturbation components in burned 
gas parallel to x it 6,2a, 6 , 2 a, ^-coordinate axes, 
respectively 

mean velocity of burned gas 
mean velocity of burned gas in turbulent 
combustion 

lateral component of velocity perturbation 
(Wi+iW 3 )/A sin <p 

amplitude coefficients of longitudinal veloc- 
ity component associated with pressure 
wave in combustible mixture 
lateral component of velocity perturbation 
(component parallel to plane of unper- 
turbed flame front and normal to u' and 
v' components) 

amplitude coefficients of lateral velocity 
component associated with pressure wave 
in combustible mixture 
position vector 

coordinate in r^y-system measured in direc- 
tion of unperturbed flame-front travel 
relative to which burned gas is at rest 
coordinate in 2 c ,y-system measured in direc- 
tion of unperturbed flame-front, travol 
relative to which combustible mixture is 
at rest 

coordinate in x u 4,23, 6-system measured 

in direction of unperturbed flame-front 
travel relative to which burned gas is at 
rest 

coordinate orthogonal to x ub and 2 3 , t and 
making angle 6 with ^-coordinate 
coordinate orthogonal to x i§ b and 22, b 
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X>,c 


*2.< 


* 3 ,e 

V 

dZ t (k,t) 


s 

a u a 2 . . . a 6 

r 


y 

A 


5i,5j . . . 5a 
«i,fa • • • eio 

r 


V 


A 

v 

/ 


coordinate in x lw c ,a^ f c -system measured 
in direction of unperturbed flame-front 
travel relative to which combustible mix- 
ture is at rest 

coordinate orthogonal to x hc and x Xc and 
making angle 6 with ^-coordinate 
coordinate orthogonal to x u c and % c 
coordinate orthogonal to x c and x b 
random amplitude vector of shear-field 
Fourier component 
coordinate orthogonal to x c and y 
grouping defined in eqs. (B18) 

-JyR t 

ratio of specific heats 

r 2 (r-l) 1 tanV+2r(T 2 +2r-l)tanV+ (r+1)* 
grouping defined in eqs. (B24) 
grouping defined in eqs. (B21) 
variable upon which pressure wave depends, 
eqs. (7) 

variable upon which pressure wave depends, 
eqs. (7) 

angle between polarization plane of incident 
shear wave and x u 0 ,xx c -plano 
coefficient defined in eq. (3e) 
k\ x„+k' s y 

flame-front displacement 
flame-front displacement velocity 



6 

P 

p' 

<T 

r 

<P 


* 

Q 

Subscripts: 

b 

c 

cr 

V 

8 

1,2,3 

i,n 

Superscripts: 

I, II 
* 

/ 


flame-front slope with respect to x 2i c -coordi- 
nate 

flame-front slope with respect to x 3: c -coordi- 
nate 

flame-front slope with respect to ^-coordinate 
static density 
static-density perturbation 
k[Ut-\-Ky, eq. (lOj) 
flame-temperature ratio, T 3 , <,/ T s . c 
spectral density tensor 
angle between wave-number vector of inci- 
dent shear wave and direction of unper- 
turbed flame-front travel, tan <p=k'Jk[ 

( ij—v ) 

scalar amplitude function defining spectral 
density tensor 

burned gas 
combustible mixture 
critical 

potential-flow velocity component 
shear-flow velocity component 
orthogonal coordinate designation 
designates order of steady-flow quantity 

designates order of fluctuation quantity 
denotes complex conjugate 
denotes fluctuation quantity except where 
otherwise specified 


APPENDIX B 


ATTENUATING- WAVE SOLUTION FOR SINGLE- WAVE INTERACTION 


The arguments of the various fluctuation quantities are 
equal at the flame front where x c =Ut and x b =(U—V)t. 
Therefore, substituting equations (10) into equations (3) 
and (5) and separately equating the real terms and the 
imaginary terms provide the following set of equations: 


sin 

P 2 +I 2 =X 2 +| EJ3 2 


(B4) 


— B 2 A sin <p— PjWj+PsZ?! 
J-i BiNi— B 3 G 2 =— BiWi+BzIli 

Ji (l~ S )-L 1 -K l N l -K l 0 1 =-K i W 1 - 

K*A sin v+K 2 B 1 +K 3 t'-(K 1 -K<)E 1 

J 2 (l- l ^-L 2 -K l N i -K l Q 2 = 

-K<W 3 + 


y (bi) 


y (B2) 


fh ^ U~—^) ^ 1 2 — to2 ( fbJiAbtJs) 

d^N i f b -jj— ^ N 3 =m 2 ( b b Ji f b J 2 


(B5) 


Pi—dtP 2 =m 2 c 1 , J 2 
dbPi—ft, jj V '^Pz=—m a c b J i 


(B6) 


— Wi-^A sin p+^Pi 

\ J 3 +L 2 +(D 1 -l)H 2 -D 1 N i -D 1 G 3 =-W 1 +± R a 




-h 


( u \ 

G% = ^ 2^2 

\TJ-V) 

- 

( U \ 

Gi=iC2ll 

> 

\u-v) 


y (B3) 


(B7) 



614 


REPORT 1277 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


/eWj— (f 0 Ri+b e R a ) 

+f c W a =^p (— b c Ri+f B R^) 


fcXi dpX. 2 — ' C c -Z?2 

^c^L+/c^s= — ^Ja CcRi 


H x — Wi— A sin ^=Al?! I 

Ho — W a =AR a J 

From equations (B5), (B6), (B8), (B9), and (BIO), 

H^=hi(h^J 2 ) J 

P i=Aj( A 4 J 1 A £ t/ 1 ) ') 

Pi = hi(hsJi j 

Hi= A £ (A t/i! i - t h^Ri) 

Wa=Ag( J 

Xi—hs{ Ao^id^AioSa 
-^ 2 = Afl(' /iio/f 1 1 

flx= A sin <pAhdJbR\~\~h&R , t) 'j 

Sa=Al2ji+Ati( — Ag-Rj+Ayiy J 


where 




wr 


n/Di+di 


A*= 


yM\f c +<%) 


n 


h^jj—h b Db A 7 ^fl—bcdc 

h^j~+dJ b h^b c f c +d c f c 


hi^Cj/Ib 

T _ 
% Q ~ 


£>! 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 


(B13) 


(B14) 


(B15) 


(B16) 


Ag = Ccd c 
Aio m c e / e 
m’mtflyO* 

From equations (B4), (B7), (B12), (B14), and (Bl5), 

ajJJj— a 3 Jj— or 4 Jj+asA sin | 

G r 2 =a a Fi+o: 1 fi s +o: 4 J 1 — 1X3J2 j 


(B17) 


where 


“■ “E& (A+^+A^,] 


a; , , 

_ K , . 

= TTFT* 


(B18) 


*“A;zV 

5= aL 335 * 

Then, from equation (B2), 

z 1 =</ 1 [i-~z- 1 (A 1 a s -^)]+j s [5: 1 (a 1 a s +« 4 )]+ 

F 2 K(AflA 8 + fl£ 2 )]+^ sin p[l£i(l— a 6 )]— -ST,/ 

£»= — Ji[^i(AiA3+a4)] + J 2 jjL — — — K,(llJh — a 3)J — 

Ei[Z , 1 (A fi A 8 +«2)]+-B2K(A 8 A 7 -«i)--K'2+ 

(Zi-ZJAJ 

The various disturbance amplitude coefficients of equa- 
tions (10) have now been obtained in terms of the coefficients 
j?i, i? 2 , Ji, J t and the parameters A sin <p and r' . From 
equations (Bl) and (B3), 


(B19) 


where 


JitiPRitzARiti — £ s 

— j 1 c/a ^1 Fi€4-j--Ka c 3 = 0 

Ji 69 -)- Jj tn~\~R 1 €g-f- R z £ 9 = £j 0 

— Bj€9-(--B2£8 = 0 J 


(B20) 


«i = 1 + B 1 — ( J5 iK 1 + B 3 ) ( At A* — “ 3 ) 

«a s (Bi-BTi-f-BaXAiAa-l-^q) 

e 3 =(BiK 1 +B 3 )(hiJi 7 —ai)+B 1 (Ki—K t )A~B3—BiK2 
£ 4 = (Bi^ti+Bs)(AjA8+“3) 

* 6 =BiK z t'~ A sin <p[(B X K X + A? 2 ) (1 — aj] 

tj = 1 — (Z?i-{--Ki)(AiA 2 <* 3 ) 

« 7 =(Z?i+^i)(A l A 3 +a4) 

€8 e=(Z? 1 +^ 1 )(A^-«i)+(2? 1 -l)A-Z- J +i+(F: 1 -X4)A 

e 9 = (Zh+liLiXAsAs+aa) 
tio^ K z t'~A sin < s(T> 1 +Z'i)(l— a 5 ) 


r (B21) 
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Ji= 


Rl ( tl*9 €469) €369) 

€3 €9 6167 


r R 1(6369 — €467) Rti&tts €367) 

WJ — 

«3«8 6467 


(B22) 


E 1= 


5j54 5i5s 


Tv 2 — 


5i5s 5^4 


S 284 5i5s 


(B23) 


where 

5 i — ei(<i«8 €4 €9) "b €3(6369 €167) €3(6467 £3*0) 

— € 3 €9 ) €4(6369 61*7) €3(6367 €363) 

^ 3 ^ € 5(6369 * 1 * 7 ) 

^ 4 3 te( *i «9 €469) 4 " €3(6369 €167) €7(6467 6469) 

Sj^6s(6i€8 6369) €7(6367 6363) 69(6969 — 6467) 

8e — ejo( 6369 6 i 67 ) 

Equations (Bll) to (B15), (B17), (B19), (B 22 ), and (B23) 
provide the formal solution of equations (Bl) to (BlO). 


V (B24) 


APPENDIX C 

ATTENUATING- WAVE SOLUTION FOR CONSTANT-PRESSURE COMBUSTION 

and 


If terms of order M 2 are retained, equations ( 1 ) provide 
the following relations for (he unperturbed-flow quantities: 


25=l-7(r-l)M 2 4- • • • 

Pc ■ 


P . U-V 
Pb ~ U 


U 1 


=r[l+^(r-l)M 2 + . . .] 

c P T s .c 1 , 1 

IP “2 + (y-1 )M 2 

2M* 


a? r[2-(y-l)(r-l)M 2 -(.y 2 -l)r(r-l)M*+ . . . ] 
With these relations, equations (4) take the form 
B l c=7rM i , Bi=2r/M 2 , B 3 =l+[l+7(r-l)]il4 J 

Dl= \[ lJ1 ir ( t “ 1 ) MS ]’ -Ex=(t- 1) [l+^ rM ! ] 

K\ =Ki= (y— 1)M 3 , AT* 

[i+2=i (t _ 1)M2 J, ^3=4l=i M 2 J 

and equations (9) may be written 
b c = b b —~k[M a , c e =c b =K, d c =m+m, d b =m+TAP ) ' 


(Cl) 


n=(Ky- p*L_iW] m* 


KC 2 ) 


If only the leading terms in powers of TVf s are retained, 
equations (B16) and (B18) provide the following relations: 


l^hi- 


yM 2 ( 1 +t 2 tan 3 <p)‘ 

khz hih 

y 

sin 2 (p 
'^yM 1 ' 

7 7 77 sin <P cos v 

fltfls /I&Iq yM* 

* 


(03) 


a ‘ = I g g [l-(r-l) tan 2 <p] 


«2= 


«S= 


r 2 sin 1 <p tan <p 
~~yM a 

t 2 tan 2 p 




yM 3 ( 1 +t 2 tan 2 <p) 


a *~yM*( 1+P tan 2 <p) 
«6=t[ 1 — (r — 1 ) tan 2 <p] J 


tan 8 1 


(C4) 


From these relations: 


7 i._. _ T tan <p l 1 _ _(t— 1 ) sin 2 tan 2 ?>— 1 ) 

hiflz~r a i y]\p* ’ a i YrVf 2 

T . . , 7 1 sin <p cos K 1 + P tan 2 ip) 

^- 03 = 0 , ^ 3 + 03 =— z— — ^ M2 — 


From equations (B21) and (C5), 


sin <p cos <p(l +T 2 tan 3 
61 «•“ 64£ «- yM^r 

(r— 1 ) sin 3 >p{t tan 2 1 ) 

€ 163 - 6369 - ^r r 


y (05) 


6367 — ejes — 6369 6167 = 


tan <p n 

yM 2 ’ e ‘ e r~ eje »— u 


( 06 ) 
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and from equations (B21), (B 24 ), and (C6): 


5i= si^s^ [( T+ i)(l+ T tan* <p)— 


2 tan* <p{i — l)r(r tan* <p — 1)] 


K t— 1 )( t tan * <P~ 1 )+ 2 (t+ 1 )t(t tan* p+ 1 )] 
5 3 =— tan ip[r ' — 2(t — l)(r tan* <p~l)A sin 90] 


5 4 = 


— (r— 1) tan <p sin*ff(T tan* ip— 1) 
fMV 

, (r+l)(r tan* 9H-I) sin* .p 

y^M*r 


y (c?) 


5 »=~lK [t' (t l)(r tan* ip~l)A sin^ 


RV-- 


Since for constant-pressure combustion R {,) = Ri/yM 2 , R w 
—RilyM 2 , and so forth, the coefficients R a) , R <3) , J m , and J (2> 
are obtained from equations (B22), (B 23 ), (C6), and (C 7 ). 
With these coefficients determined, the remainder are 
obtained from equations (B 7 ), (Bll) to (B 14 ), (B 17 ), and 
(B 19 ). The amplitude coefficients are 


R « 


Ni- 


-^(r — l)*(r tan* tp — 1)*A sin <p — [t*(t*-[-2t — 1) tan* ¥>+t(t*+1)] — 

= J (J, = — — — -^(t* — l)(r* tan 4 cp—l)A sin <p-|- [ t 3 (-3 — 1) tan 4 <p~ 2 -r 3 tan* <p— t(t+1)] ^7^- 

ja) = _l. tan* sin H-t*[t(t-1)* tan 4 ^+(t*+2t-1) tan* <p+ 2] ^ j- 

L\= j L 2 — 0 

T 

A(l+r*tan a y ) (( T— 1 )( T fcaa * 1 )[(t+1)(t tan* p+1 )— t(t— 1) tan* <p(r tan* *>— 1)]A sin <p— 
t-^ r* tan* v>[t(t — l) 2 tan 4 ^>+t(t+ 1) tan* 9H*4]+(r+l)^-~^ 

{ (t tan2 *- 1 )[ t(t+ 1 ) (t tani H-l)+(r-l)(r tan* ip-l)\A sin ^+ 
r*[r(2i — 1) tan 4 <p — (r — 1) tan* ip — 1] — 

Ql= A( r +r* r tan* <p ) K T +1)*( T tjmS 1) 2 H-t (t — l) 2 tan* <p(r tan* <p— 1)*] J. sin <p— 

t ( i — 1) tan* ¥?[t(t*-(-2t — 1) tan* ^+(t*-!-1)] — 

G 3 = qj^tan*^ tan* <p— l)(r tan* <p+l)A sin <p+[r 3 (r— 1) tan 4 <p— 2r 3 tan* <p— t(t+1)] 7^7 J- 

Fi=— iV 2 , P 2 =Ni 

t Gi r —G 3 
A= — 7— — ’ 


W i 


r tan <p t tan ip 

\= — i -^(r — 1 )( T tan* ip — 1 )[t(t — 1 ) tan* (p-bOr-f-lH-A sin <p — [t*( 3 t — 1 ) tan* <p+t(t-(- 1 )] 

— — P tan jg -^2(r tan* 97—1)44 sin p+t[t tan* p+1] — 

Xi=&2 = Wi, X 3 — — Wi 


W 3 -- 


Sj=— -^2(tt)-1)(t tan* ip+l)^ sin ipA[t( 3 t — 1) tan* yd-fr-l-l)] — 
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For constant-pressure combustion, Sc^Zb—k'JJt+kiy, 
y e =k' 2 (x e -Ut), Vb =k' 2 [-x b +(U-V)t], y=j c ' 1 x e +k' 2 y,t= 1 ^+ 

T 

k 2 y, and <x=k[Ut-{-k 2 y. At the flame front where x c = Ut 
and x b =(U—V)t, y c =Vb=0, the attenuation factors e~** 
and are 1, and the arguments of the disturbance waves 
have the form (k[ Ut+kUj) . 

The procedure of reference 23 is utilized in reference 22 to 
treat a similar interaction problem for consideration of 
flame-front stability. It is assumed that there are first- 
order perturbations in the laminar flame speed U (in the 
present notation). Reference 22 presents results only for 
the special case where the plane of polarization of the inci- 
dent shear wave is in the plane x Uc —x^ e (8=0°) and the 
wave-number vector is parallel to the unperturbed 'flame 
front (<p=90°). The results of the present analysis were 
compared with those of reference 22 for the case of an 
absence of first-order perturbations in the laminar flame 
speed and for the special case of 6=0° and ^=90° without 
heat-release perturbations. Although agreement as to sig n 
and magnitude is obtained for the shear-field velocity com- 
ponents, apparently differences in sign occur for the potential- 
field velocity-component amplitudes as indicated in the 
following table: 



Present 

analysis 

Analysis of 
re£ 22 

V 

-A 

~A 

XJ 

IA 

-IA 

ok 

u 

A 

-A 

f 

U 

iA 

iA 


Consideration of conservation of momentum across the 
flame front indicates that v' v ,dU and v' v , b jU should be of the 
same sign. The present analysis is in agreement with this 
consideration. 
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